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Abstract 

Time correlated fluctuations interacting with a spatial asymmetry po- 
tential are sufficient conditions to give rise to transport of Brownian par- 
ticles. The transfer of information coming from the nonequilibrium bath, 
viewed as a source of negentropy, give rise to the correlated noise. The 
algorithmic complexity of an object provides a means of quantitating its 
information contents. The Kolmogorov information entropy or algorith- 
mic complexity is investigated in order to quantitate the transfer of in- 
formation that occurs in computational models showing noise induced 
transport. The complexity is measured in terms of the average number 
of bits per time unit necessary to specify the sequence generated by the 
system. 
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Recently a lot of attention has been paid to simple stochastic models in which 
some of the energy in a noncquilibrium bath is used to obtain a net transfer of 
information. The best known examples of this process are thermal ratchets or 
correlation ratchets where a nonzero net drift speed may be obtained from time 
correlated fluctuations interacting with broken symmetry structures. Recently 
the study of a wide variety of ratchets has received much attention due to their 
general interest and possible technological applications in modeling molecular 
motors and other nanoscale and mesoscale systems. For a general review of 
thermal ratchets the reader is referred to Ref. jjj . 

As was demonstrated by M. Millonas |2j, the net current is obtained by 
means of a source or sink of negentropy (physical information) that allows the 
engine to operate. Numerous existing studies have addressed the question of 
the efficiency in ratchets || but, despite this widespread interest none of them 
did it from the perspective of information theory and, up to our knowledge, a 
quantitative estimate of the amount of negentropy was never obtained. A possi- 
ble reason that no effort has been done to study the entropy transfer on thermal 
ratchets may be that they constitute far from equilibrium, non-ergodic systems 
and usual entropy calculations rely on probabilistic ensemble concepts that re- 
quire energy minimization procedures associated with equilibrium arguments. 
Algorithmic complexity was a concept introduced by Kolmogorov to avoid the 
probabilistic interpretation of Shannon that has been recently used to obtain 
the information content of far from equilibrium systems such as proteins and 
fractal growth processes [Q. This characteristic makes algorithmic complexity 
specially useful to obtain the variation of entropy or physical information of far 
from equilibrium systems such as thermal ratchets. 

Following the idea developed by Crisanti et al || and Paladin et al [|| that 
applied algorithmic complexity to characterize the chaoticity of dynamical sys- 
tems with noise, we study the algorithmic complexity of thermal ratchet motion. 

The algorithmic complexity of an object is broadly defined as the length in 
bits of the shortest description for that object In other words, complexity 
can be characterized as the length of the shortest possible definition of the 
sequence itself. 

In the framework of information theory, if one wants to transmit the se- 
quence, through a noisy channel, until it deviates from the deterministic one a 
tolerance threshold A, one can use the following strategy: 

Specify the initial condition x(l) with precision So using a number of bits 
n = log2{A/5o) which permits to arrive up to a time n where the error equals 
A. Then specify again a new initial condition x(ti + 1) with a precision Sq 
and arrive up to a time T2, and so on N times. The number of bits necessary 
to specify a sequence with a tolerance A up to T max = ^t=i T i i s — Nn and 
the mean information generated per time step is ~ Nn/T max = K a /ln2 bits. 
Following the idea of algorithmic complexity as the length in bits of the mean 
information, K a is the complexity associated to stochastic motion. Thus, K a 
can be obtained as M: 
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where 

1 - 

t= lim — Vt,. (2) 

i=l 

The model for ratchets that we have used is based on the model proposed 
by Astumian and Bier 0. Two of the authors, recently used this model to 
study the approach to steady state in ratchets || and to study ratchets with 
finite inertia^]. It consists of an asymmetric piecewise linear potential where 
the barrier height fluctuates between two states. Astumian and Bier showed 
that the fluctuations of the barrier height can produce a net flow even with the 
net macroscopic force that is zero at all times. The general Langevin equation 
describing the ratchets that we will consider is 

dx 

— = -u'(x,t) + y/2k B Tw(t), (3) 

where u'(x, t) — ^jff^-, T is the temperature and k B is the Boltzmann constant. 
The potential u{x, t) switches on and off the v{x) asymmetric potential, and w(t) 
represents the Gaussian noise term with delta function correlation {w(t)w(t')) = 
5{t-t'). 

The fluctuations cause the potential u(x,t) to switch between u+(x) = v(x) 
and U-(x) — with a dichotomous Markov process governed by the master 
equation 

d f P +(t)\ _J~l 1 \ (P+ 



dt \P-(t) 
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where 7 represents the frequency of the flipping between and U—(x). 

In our numerical simulations we have followed the approach of Elston and 
Doering . Details of the simulation procedure can be found in references Q 
and§. 

In order to obtain the algorithmic complexity associated to ratchet motion, 
two different walkers arc allowed to move under the ratchet potential with initial 
conditions separated by So, x(l) and x'(l) = x(l) + Sq. Both motions continue 
until the difference equals A at time t\ . Then begin again with initial conditions 
x(ti + 1) and x'(ti + 1) = x(ti + 1) + So, and so on N times. 

Simulations were carried out on ratchets with an asymmetric barrier |?], || 
switched on and off with different flipping rates. The asymmetric parameter 
used was a = 10/11 and the barrier height was E = 10. This is a very im- 
portant kind of ratchet because the experimental studies of ratchets done by 
Rousselet et al |llj and by Faucheux et al Jl^] were carried out with an asym- 
metric potential switched on and off periodically and Ajdari and Prost [|l3) used 
the on-off barrier fluctuation as a method for dielectrophosphoretic separation. 

In the limits of low and high flipping rates the current tends to zero and 
the system behaves essentially as a random walk. The relative complexity for a 
given flipping rate is defined as: 

K - K RW 
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Figure 1: Comparison between absolute value of the current (solid 
line) and relative complexity (circles) of an on-off ratchet vs 
logarithm of flipping rate. 



where K is the complexity at the given flipping rate and Krw is the complexity 
for high flipping rate limit. The relative complexity and the net current as 
functions of the logarithm of the flipping rate 7 are shown in Fig. 1. In this 
figure it can be appreciated that the relative complexity tends to zero in the 
limits of low and high flipping rates, where the current tends to zero. On the 
other hand a dramatic increase of the relative complexity appears for flipping 
rates corresponding to the highest currents. 

The proportionality between K and J can be explained using the well known 
statistical mechanical relationship derived by Zurek |T3 



where S is the physical entropy, K is the algorithmic complexity and / is the 
Shannon information entropy. Zurek's relationship says that the physical or 
thermodynamic entropy of a system is composed of two parts, that determined 
from the known information of the system K and that determined from the 
unknown or probabilistic information /. When the microstates of a system are 
unknown K = and the entire entropy is due to the Shannon information, i. e., 



S = K + 1, 



(6) 
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Figure 2: Log-log plot of the relative complexity vs the error A 
for high and low current flipping rates. 



S = I. As observations are made on the system the information content shifts 
from / to K. 

When there is a net current in the ratchet, the amount of known information 
about the system increases with respect to the case of no net current. The exis- 
tence of a net current implies there is some more information known about the 
system. This amount of extra information comes from the negentropy supplied 
by the external nonequilibrium bath. As a consequence, the maximum amount 
of extra information or complexity is obtained when the current is maximum. 

Another way to characterize the behavior of complexity in on-off ratchets is 
the scaling of the complexity with the error A. We studied the scaling behavior 
of complexity for the random walk limit corresponding to a high flipping rate 
and for the highest current flipping rate. The scaling behavior for both cases 
is shown in Fig. 2 where log(K) vs log(A) is plotted. In the random walk 
limit K <~ A~ 2 according to the standard diffusion law. For the highest current 
flipping rate, due to the existence of a net current, the complexity scales as 
K ~ A -1 . 

The simulations were carried out using a rather high value of A ~ 10 in order 
that the interaction between the ratchet potential and the on-off fluctuation may 
take effect. Results were averaged over 500 independent runs. 
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In conclusion, the net transfer of information from a source of negentropy 
that allows the net drift of particles in thermal ratchets may be obtained through 
the algorithmic complexity. The negentropy transference results in an increase 
of complexity observed when high ratchet currents take place (see Fig. 1). Ad- 
ditionally a different scaling behavior of complexity at high currents compared 
with the no current case, i.e. random walk limit, is observed as can be seen in 
Fig. 2. As far as we know, the algorithmic complexity lets us obtain, for the 
first time, an estimate of the net transfer of information associated with ratchet 
motion. 
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